Theory of pump-probe experiments of metallic metamaterials coupled to the gain medium 
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We establish a new approach for pump-probe simulations of metallic metamaterials coupled to the gain ma- 
terials. It is of vital importance to understand the mechanism of the coupling of metamaterials with the gain 
medium. Using a four-level gain system, we have studied light amplification of arrays of metallic split-ring 
resonators (SRRs) with a gain layer underneath. We find that that the differential transmittance AT /T can be 
negative for SRRs on the top of the gain substrate, which is not expected, and AT/T is positive for the gain 
substrate alone. These simulations agree with pump-probe experiments and can help to design new experiments 
to compensate the losses of metamaterials. 

PACS numbers: 42.25.-p, 78.20.Ci, 41.20.Jb 
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The field of metamaterials has been spectacular experi- 
ments progress in recent years |1]J-|3|. The mostly metama- 
terial are metal-based nano-structure and eventually suffering 
from the conductor losses at optical frequencies, which are 
stills orders of magnitude too large for the realistic applica- 
tions. In addition, metamaterial losses become an increas- 
ingly important issue when moving from multiple metal-based 
metamaterial layers to the bulk case p| . Thus, the need for re- 
ducing or even compensating of the losses is a key challenge 
for metamaterial technologies. One promising way of over- 
coming the losses is based on introducing the gain material to 
the metamaterial. The idea of combination of a metamaterial 
with an optical gain material has been investigated by sev- 
eral theoretical |J4H7|] and experimental studies ll8l4l2ll. From 
the experiments point of view, the realistic gain can be ex- 
perimentally realized with fluorescent dye s llsll, qu antum dots 
ISl llOD or semiconductor quantum wells lllli 11211 . All these 
loss-compensation are mainly attributed to the coupling be- 
tween metamaterial and the gain medium. Without sufficient 
coupling, no loss-compensation can happen, nor can the trans- 
mitted signal be amplified. Therefore, it is of vital importance 
to understand the mechanism of the coupling between meta- 
material and the gain medium. In addition, these ideas can be 
used in plasmonics to incorporate gain lll3l 11411 to obtain new 
nano-plasmonic lasers il5lll6ll . 

In this Letter, we present a systematic theoretical model 
for pump-probe experiments of metallic metamaterials cou- 
pled with the gain material, described by a generic four-level 
atomic system. We describe the dynamical processes in meta- 
materials with gain, and increasing the gain changes the meta- 
material sproperties and we need to have self-consistent calcu- 
lations (|4|-|6[] to reach a steady state. The pump-probe results 
affecting the time dependence of the population inversion and 
the electric field enhancement that increases the effective gain. 
We observe differential transmittance signals from the cou- 
pled system that are larger than for the bare gain. Further- 
more, we observe a more rapid temporal decay of the differ- 
ential transmittance signal for the coupled system compared 
to the bare gain. Both effects indicate substantial local-field- 



enhancement effects, which increase the effective metamate- 
rial gain beyond the bare gain, leading to a significant reduc- 
tion of the metamaterial's losses. 

We model the dispersive Lorentz active medium using a 
generic four-level atomic system. The population density in 
each level is given by Ni (i=0,l,2,3). The time-dependent 
Maxwell's equations for isotropic media are given by V x 
E(r,i) = -dB{r,t)/dt and V x H(r,i) = dT){r,t)/dt, 
where B(r,t) = /x/XoH(r,i), D(r, t) = esoE{r,t) +P{r,t) 
and P(r, t) is the dispersive electric polarization density that 
corresponds to the transitions between two atomic levels, A^i 
and N2. The vectors P introduces gain in Maxwell's equa- 
tions and its time evolution can be shown to follow that of a 
homogeneously broadened Lorentzian oscillator driven by the 
coupling between the population inversion and external elec- 
tric field iIitII . Thus, P obeys the equation of motion 
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-cj2p(r,t) = aaAN{r,t)E{r,t) 



where Ta stands for the linewidth of the atomic transitions at 
cJq, and accounts for both the nonradiative energy decay rate, 
as well as dephasing processes that arise from incoherently 
driven polarizations. In the following simulations, this value 
is equal to 27r x 20 • 10^^ rad/s. o-q is the coupling strength 
of P to the external electric field and its value is taken to be 
10-4 CVkg. The factor AN{r, t) = Ni{r, t)~N2{r, t) is the 
population inversion between level 2 and level 1 that drives 
the polarization P. In order to do pump-probe experiments 
numerically we first pump the gain material with a short, in- 
tense Gaussian pump pulse. After a suitable time delay we 
probe the structure with a weak probe pulse (see Fig. [T). In 
our model, an external mechanism pumps electrons from the 
ground state level Nq to the third level A'3 using a gaussian 
pumping Pg{t), which is proportional to the pumping inten- 
sity in the experiments. After a short lifetime T32 electrons 
transfer non-radiative into metastable second level N-z- The 
second level (iV2) and the first level (iVi) are called the upper 
and lower lasing levels. Electrons can be transferred from the 
upper to the lower lasing level by spontaneous and stimulated 



pump pulse 



population inversion 




FIG. 1: (Color online) Schematic illustration of pump-probe experi- 
ments. 



emission. At last, electrons transfer quickly and non-radiative 
from the first level (iVi) to the ground state level (iVo). The 
lifetimes and energies of the upper and lower lasing levels are 
T21, £'2 and Tio, El, respectively. The center frequency of 
the radiation is Wa = {E2 — Ei)/h which is a controlled 
variable chosen according to the pump-probe experiments. 
The parameters T32, T21, and tio are chosen to be 0.05ps, 
80ps, and 0.05ps, respectively. The initial electron density. 



Noir,t = 0) = 5.0 X 10 
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(i=l,2,3). Thus, the atomic population densities obey the fol- 
lowing rate equations: 
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where Gaussian pump Pgit) = Po x e 



-f 



Tio 



with Pq = 



3 X lO^s" 



ti) 



6ps 111811 . and Tp = 0.15 ps. 



In order to solve the response of the active materials in 
the electro mag netic fields numerically, the FDTD technique 
is utiUzed 11911 . using an approach similar to the one outlined 
inQ. 

The object of our studies is to present pump-probe sim- 
ulations on arrays of silver SRRs coupled to single quan- 
tum wells illi I12I1 . The structure considered is a U-shape 
SRRs fabricated on a gain-GaAs substrate with a square pe- 
riodicity of p = 250nm (see Fig. |2ja)). The SRRs is made 
of silver with its permittivity modeled by a Drude responds; 
e{uj) = 1 - cj^/(w^ + iujj), with ujp = 1.37 x 10^'^ rad/s and 
7 = 2.73 x 10^'^ rad/s. The incident wave propagates perpen- 
dicular to the SRRs plane and has the electric field polariza- 
tion parallel to the gap (see Fig.|2a)). The corresponding geo- 
metrical parameters are a = 150nm, hd = 40nm, hg ~ 20nm, 
hs = 30nm, w — 50nm, and h = 75nm. Fig. |2b) shows 
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FIG. 2: (Color online) (a) Schematic of the unit cell for the silver- 
based SRRs structure (yellow) with the electric field polarization par- 
allel to the gap. The dielectric constants e for gain (red) and GaAs 
(light blue) are 9.0 and 11.0, respectively, (b) Calculated spectra 
for transmittance T (black), reflectance R (red), and absorptance 
A (blue) for the structure shown in Fig. lUa). The insert shows 
the profile of the probe pulse with a center frequency of 175THz 
(FWHM = 2THz). 



the calculated spectrum (without pump) of transmittance T, 
reflectance R, and absorptance A for the structure shown in 
Fig. |2a). The resonant frequency is around 175THz, and 
we refer to the resonant frequency according to the dip of 
the transmittance. In our analysis, we first pump the active 
structure (see Fig.|2ja)) with a short intensive Gaussian pump 
pulse, Pg{t), (see Fig. [3] top panel). After a suitable time de- 
lay (i.e. the pump-probe delay), we probe the structure with 
a weak Gaussian probe pulse with a center frequency close to 
the SRRs resonance frequency of 175THz. Typical examples 
for the spatial distribution of electric field an gain are shown in 
fl2Ill . The incident electric field amplitude of the probe pulse 
is lOV/m, which is well inside the linear response regime. 
Then, we can Fourier transform the time-dependent transmit- 
ted electric field and divide by the Fourier transform of the in- 
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FIG. 3: (Color online) Schematic of the numerical pump-probe ex- 
periments for the case on resonance. From the top to the bottom, 
each row is corresponding to the pump pulse, population inversion, 
incident signal (with time delays 5ps, 45ps and 135ps), transmit- 
ted signal, and differential transmittance AT/T. It should be men- 
tioned here the incident frequency of the probe pulse is 175THz with 
FWHM of 2THz and is equal to the SRRs resonance frequency. 



cident probe pulse to obtain the spectral transmittance of the 
system as seen by the probe pulse. Additionally, we obtain the 
total pulse transmittance by dividing the energy in the trans- 
mitted pulse by the energy in the incident pulse, integrated in 
time domain. We define the differential transmittance, AT/T, 
by taking the difference of the measured total plus transmit- 
tance with pumping the active structure minus the same with- 
out pumping and dividing it by the total plus transmittance 
without pumping. This differential transmittance is a function 
of the pump-probe delay. The bottom panel in Fig. [3]gives a 
differential transmittance AT/T which is negative, this result 
was not expected, and we need to understand this behavior, 
which agrees with the experiments lllUll2 l. 

Fig. m gives an overview over the results obtained for the 
case of the SRRs on resonance, i.e., w^ = 27r x 175 x 10^'^ 
rad/s. Data for the structure in Fig. |2ja) (left column in Fig.|4]i 
and for the bare gain case (right column in Fig. H) without the 
SRRs on top is shown. For parallel polarization, the light does 
couple to the fundamental SRRs resonance, for perpendicular 
polarization it does not. The probe center frequency decreases 
from top (179THz) to bottom (169THz). Note that the width 
of the probe spectrum is 2THz (see the insert of Fig. EJa)). 
Hence, the data have been taken with 2-THz spectral separa- 
tion. Inspection of the left column shows a rather different 
behavior for the SRRs with gain compared to the bare gain 
case. While the bare gain always delivers positive AT/T sig- 



nals below +0.16% (right column) over the whole probe spec- 
trum. The sign and magnitude of the signals change for the 
case SRRs with gain. Under some conditions, AT/T reaches 
values as negative as —8.50% around f probe = 175THz. Ad- 
ditionally, we may also get positive AT/T at the very edges 
of the probe range (see left column in Fig. |4|i. If we turn to 
the case of perpendicular polarization case, no distinct change 
between the pump-probe results on the SRRs (not shown in 
Fig. IHi and the bare gain (right column in Fig. |4]i, neither in 
the magnitude nor in the dynamics of the AT/T, can be de- 
tected. 

We argue the distinct behavior can be attributed to the 
strong coupling between the resonances of the SRRs and the 
gain medium. The negative AT/T are not as we expected 
at first glance: the pump lifts electrons from ground state 
to an excited state so that the absorption of the probe pulse 
is reduced, leading to an increase of transmission. This is 
not the whole story. The reason lies in the fact that with 
pump we not only affect the absorption, but disturb the re- 
flection of the structure, resulting in the mismatching of the 
impedance. Furthermore, we observed either increasing or 
decreasing tendency for the case of on resonance as shown 
in Fig. |4] All those behaviors can be explained by the com- 
peting of the weak gain resonance and the impedance mis- 
matching between pump and without pump cases. We will 
explore the underlying mechanism below. Fig. |5] shows the 
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FIG. 4: (Color online) Time domain numerical pump-probe experi- 
ments results for the SRRs that is nearly on-resonant with the gain 
material. The left column conesponds to the parallel probe polariza- 
tion with respect to the gap of the SRRs, the right column is the case 
for bare gain material, i.e., without SRRs on the top of the substrate. 
The width of the probe signal is 2THz with decreasing in the probe 
center frequency from 179THz for the top panel to 169THz for the 
bottom panel. 

results for the difference in absorptance (AA), difference 
in reflectance (AR), their sum (A^ + AR), and the dif- 
ference in transmittance (AT = —{AA + AR)) between 
pump (Po = 3 X 10^s~^) and no-pump using a wide probe 
(FWHM = 54THz) pulse with a fixed pump-probe delay 
of 5ps. As expected, we may observe a positive differential 
transmittance, AT/T > 0, when we pump the gain, A^ < 0, 



and if AR (impedance match) remains unchanged. 

The results of Fig. |5] are obtained for pump-probe exper- 
iments with the probe frequency equal to the resonance fre- 
quency of the SRRs (175THz) at a pump-probe delay of 5ps; 
Results for longer pump-probe delays are shown in supple- 
mentaty material 12111 . Notice that AR is positive, A A is neg- 
ative, and AT is also negative very close the the resonance 
frequency. If the probe center frequency moves away from 
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FIG. 5: (Color online) Frequency domain numerical pump-probe 
experiments results for the on-resonance case. Simulations results 
for the differences in transmittance (AT), reflectance (A_R), and ab- 
sorptance (A^) versus frequency. 



the SRRs resonance frequency, the negative AT /T decreases 
in magnitude and finally AT /T becomes positive. These re- 
sults are shown in Fig. |6] and agree with experiments lllllll2 l. 
If we can increase the magnitude of the Gaussian pump pulse 



Pg{t) to 5 X lO^^^s"^ and we repeat the pump-probe exper- 
iments, AT jT « —100% at resonance frequency, 175THz. 
If we increase the pump amplitude further to lO^^s"^ we can 
compensate the losses. However, such pump intensities are 
unrealistic experimentally Bill . In conclusion, we have intro- 
duced a new approach for pump-probe simulations of metallic 
metamaterials coupled to gain materials. We study the cou- 
pling between the U-shaped SRRs and the gain material de- 
scribed by a 4-level gain model. Using pump-probe simula- 
tions we find a distinct behavior for the differential transmit- 
tance AT jT of the probe pulse with and without SRRs both 
in magnitude and sign (negative, unexpected, and/or positive). 
Our new approach has verified that the coupling between the 
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FIG. 6: (Color online) The transmittance T (without pump, solid 
line) and the on-resonance differential transmittance AT/T results 
(vector arrow). The direction and the length of the arrow stand for 
the sign and the amplitude of AT/T, respectively. The squares from 
Pi to Pe correspond to the frequency of probe pulse ranging from 
169THz to 179THz with uniform step of 2THz. 



metamaterial resonance and the gain medium is dominated by 
near-field interactions. Our model can be used to design new 
pump-probe experiments to compensate the losses of meta- 
materials. 
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